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Abstract 

A search for — > P e oscillations was conducted by the Liquid Scintillator Neutrino Detector 
at the Los Alamos Neutron Science Center using from fi + decay at rest. A total excess of 
87.9 ± 22.4 ± 6.0 events consistent with u e p — > e + n scattering was observed above the expected 
background. This excess corresponds to an oscillation probability of (0.264 ± 0.067 ± 0.045)%, 
which is consistent with an earlier analysis. In conjunction with other known limits on neutrino 
oscillations, the LSND data suggest that neutrino oscillations occur in the 0.2 — 10 eV 2 /c 4 Am 2 
range, indicating a neutrino mass greater than 0.4 eV/c 2 . 
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I. INTRODUCTION 



A neutrino produced in a weak decay is always from a specific family, v e , or is T , that 
is directly associated with the charged lepton accompanying the decay. When this neutrino 
is detected in a charged- current reaction, it manifests its identity by transforming into the 
anti-particle of the charged lepton that accompanied its creation. Lepton family number is 
then conserved. However, the result is different if the neutrino changes from one family to 
another. For example, if a changes to a z/ e , then a ji + is made at the neutrino's creation 
and an e~ created at its demise, in clear violation of lepton family number. Such neutrino 
oscillations are viewed as possible, or even likely, as the flavor eigenstates (z/ e , z/ r ) need 
not be neutrino mass eigenstates (i>i, z/ 2 , z/3). If the neutrino flavor eigenstates are a linear 
combination of the mass eigenstates, the neutrino flavor must change with time because the 
phases of the mass eigenstates evolve at different rates. In the case of two flavor eigenstates 
(z/ a , u b ), the probability that u a will turn into v h is given by 



where 9 is the mixing angle between the mass eigenstates v\ and is 2 , Am 2 is the difference 
in neutrino eigenstate masses squared, m\ — m 2 ., in eV 2 /c 4 , L u is the distance traveled by 
the neutrino in meters, and E v the neutrino energy in MeV. 

A search for neutrino oscillations requires knowledge of the neutrino source, both with 
respect to the flavor composition and energy spectrum of the source. There are two types of 
searches. The first seeks to observe a reduction in the expected number of detected neutri- 
nos of a specific flavor. Characterizing the reduction as P(aa) = 1 — P(ab), it can then be 
explained in terms of neutrino oscillations. Such searches are termed disappearance experi- 
ments. The second method looks for a greater than expected number of events ascribed to a 
neutrino flavor that is either absent or very weakly produced at the neutrino source. These 
searches are referred to as appearance measurements. The results of the search reported in 
this paper are of the latter kind. It reports an excess of events ascribed to electron antineu- 
trinos that is approximately five times greater than the number of such events believed to 
be created at the neutrino source. 

Neutrinos are assumed to be massless in the Standard Model, so the observation of 
neutrino oscillations would require an extension of the current version. In addition, as there 
are ~ 100 neutrinos per cm 3 of each neutrino family left over from the initial expansion 




(1.1) 
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of the universe, neutrino masses of even a few eV/c 2 would have a significant effect on the 
evolving structure of the universe. 

The source of neutrinos for the measurement in this report is the interaction of the intense 
(~ 1 mA) 798 MeV proton beam at the Los Alamos Neutron Science Center (LANSCE), 
which produces a large number of pions, mostly 7r + . The 7r~ are mainly absorbed and 
only a small fraction decay to fi~ , which in turn are largely captured. Thus, the resulting 
neutrino source is dominantly due to 7r + — > /i + z/ M and fi + — > e + u e i? fl decays, most of which 
decay at rest (DAR). Such a source has a paucity of z/ e , and so measurement of the reaction 
u e p — > e + n, which has a large and well known cross section, provides a sensitive way to search 
for — > v e oscillations. Such events are identified by detection of both the e + and the 2.2 
MeV 7 from the reaction np — > d'j. In addition, the v e flux from tt + and fi + decay-in-flight 
(DIF) is very small, which allows a search for —> v e oscillations via the measurement of 
electrons above the Michel electron endpoint from the reaction v e C — > e~N. 

The Liquid Scintillation Neutrino Detector (LSND) experiment took data over six calen- 
dar years (1993-1998). During this period the LANSCE accelerator operated for 17 months, 
delivering 28,896 C of protons on the production target. Using partial samples of the re- 
sulting data, evidence for neutrino oscillations has been published previously Jl], 0, This 
report presents the final results on oscillations using all the data, combining the — > v e 
and — > v e searches into a single analysis with common selection criteria, and employing 
a new event reconstruction that greatly improves the spatial resolution. An excess of events 
consistent with neutrino oscillations is observed which requires that at least one neutrino 
have a mass greater than 0.4 eV/c 2 . 

Neutrino oscillations have also been employed to explain the observed deficit of solar 
neutrinos |4[] and the atmospheric neutrino anomaly by v e and z/ M disappearance, respec- 
tively. The Super Kamiokande atmospheric results favor — > v T and provide compelling 
evidence for neutrino oscillations. It is difficult to explain the solar neutrino deficit, the at- 
mospheric neutrino anomaly, and the LSND excess of events with only three flavors of 
neutrinos, so that a fourth, sterile neutrino has been proposed to explain all of the data . 
Neutrino oscillations between active and sterile neutrinos could have a significant effect on 
the R process in type II supernovae |J . 
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II. NEUTRINO BEAM, DETECTOR, AND DATA COLLECTION 



A. Proton Beam and Targets 

The LSND experiment || was designed to search for — > u e oscillations from fi + DAR 
with high sensitivity. A layout of the detector and beam stop is shown in Fig. |l|. The 
LANSCE accelerator is an intense source of low energy neutrinos produced with a proton 
current of 1 mA at 798 MeV kinetic energy. For the 1993-1995 running period the production 
target consisted of a 30-cm long water target (20-cm in 1993) followed by a water-cooled Cu 
beam dump, while for the 1996-1998 running period the production target was reconfigured 
with the water target replaced by a close-packed, high-Z target. The muon DAR neutrino 
flux with the latter configuration was only 2/3 of the neutrino flux with the original water 
target, while the pion DIF neutrino flux was reduced to 1/2 of the original flux. The resulting 
DAR neutrino fluxes are well understood because almost all detectable neutrinos arise from 
7r + or fi + decay; tt~ and fi~ that stop are readily captured in the Fe of the shielding and 



Cu of the beam stop ||10|| . The production of kaons or heavier mesons is negligible at these 
proton energies. The v e flux is calculated to be only ~ 8 x 10~ 4 as large as the flux in the 
20 < E v < 52.8 MeV energy range, so that the observation of a v e event rate significantly 
above the calculated background would be evidence for — > v e oscillations. 

For the first three years of data taking, thin carbon targets were in place in positions Al 
and A2 at the experimental area of the LANSCE accelerator, but dominant pion production 
occurred at the A6 beam stop, which accounted for ~ 98% of the DAR neutrino flux and 
~ 95% of the DIF neutrino flux. The Al, A2, and A6 targets were approximately 135 m, 
110 m, and 30 m, respectively, from the center of the LSND detector. A6 was essentially the 
only source of neutrinos for the last three years of data taking. Note that in each case there 
was a small open space downstream of the primary targets where a few percent of the pions 
decay in flight, producing up to an energy of 300 MeV. The neutrino flux was calculated 
by a program Jl0| using particle production data for thin targets taken at a number of proton 
energies and extrapolated to the actual geometry represented. Fig. shows the layout of 
the A6 beam stop as it was configured for the 1993-1995 data taking. Table | shows the 
proton beam statistics for each of the six years of running from 1993 through 1998. 
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B. Neutrino Sources 



Neutrinos arise from both pion and muon decays. The pion decay modes are ix + — > 
7r + — > e + z/ e , 7r~ — > and 7r~ — > e~i> e . The 7r + decay occurs both with the pion at 

rest (97%) and in flight (3%). The ir~ , however, only decays in flight as they are totally 
absorbed on nuclei when they stop. Helium represents an anomalous case in which ir~ decay 
occurs occasionally, but this effect is negligible in other nuclei ||11|| . Muon decay modes are 
fi + — > e + i> e u ll and pT — > e~z/ e i/ M . Almost all fi + stop before decaying and produce a normal 
Michel spectrum for u e and v^. The yT are produced following 7r~ DIF and either decay 
in orbit or are absorbed in a nucleus through fi~N — > ^X, where .E^ < 90 MeV. The 
absorption rates are taken from [jllj and are shown in Table |J. Each of these production 
processes has been included in the flux calculation described below. 



Production Monte Carlo 



The Production Monte Carlo [10] simulates the decays of pions and muons for each of 
the decay and absorption reactions described above and for each of the configurations listed 
in Table |I[ Pion production data using a number of different proton energies were input, as 
well as information on the target materials. The particles were tracked through the specified 
materials and geometries. For each configuration, the flux and energy spectrum of neutrinos 
from each decay channel were obtained for 25 different positions within the detector. For 
DAR neutrinos the flux is isotropic. The accumulated charge of beam protons was used to 
obtain the number of protons on target, and for each year of running the resulting fluxes 
and spectra from all configurations were added together, weighted by the accumulated beam 
charges. The program gives fluxes in terms of the number of neutrinos traversing the detector 
region per proton on target per unit of area. 



D. Neutrino Fluxes 

Fig. |3] shows the neutrino energy spectra from the largest DAR sources. The flux 
from /i + DAR provides the neutrinos for the z/ M — > v e oscillation analysis. The v e flux from 
fi + DAR provides events used to verify the DAR neutrino fluxes, as discussed later in this 
paper. The v e flux from pT DAR is a background to the oscillation signal with an energy 
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spectrum similar to that of v e from /i + decay. 

Fig. ^ shows the neutrino energy spectra from various DIF sources averaged over the 
detector. The flux from ir + DIF provides neutrinos for the — > v e oscillation analysis. 
The u e flux from fi + and ir + DIF is a background for the DIF oscillation analysis. The fi + 
DIF flux is suppressed due to the long muon lifetime, while the 7r + DIF flux is suppressed 
due to the small tt + — > e + i/ e branching ratio of 1.2 x 10 -4 . 

Calculations of fi + DAR fluxes are uncertain at the 7% level, while ir^ 1 DIF fluxes and 
fi~ DAR fluxes are uncertain to 15% M. Neutrino fluxes for different years are shown in 



Table III. 



E. Detector 

The LSND detector |J consisted of an approximately cylindrical tank 8.3 m long by 5.7 m 
in diameter. The center of the detector was 30 m from the A6 neutrino source. On the inside 
surface of the tank, 1220 8-inch Hamamatsu phototubes (PMTs) covered 25% of the area 
with photocathode. The tank was filled with 167 metric tons of liquid scintillator consisting 
of mineral oil and 0.031 g/1 of b-PBD. This low scintillator concentration allows the detection 
of both Cerenkov light and scintillation light and yields an attenuation length of more than 
20 m for wavelengths greater than 400 nm ||13|| . A typical 45 MeV electron created in the 
detector produced a total of ~ 1500 photoelectrons, of which ~ 280 photoelectrons were in 
the Cerenkov cone. PMT time and pulse-height signals were used to reconstruct the track 
with an average RMS position resolution of ~ 14 cm, an angular resolution of ~ 12 degrees, 
and an energy resolution of ~ 7% at the Michel endpoint of 52.8 MeV. The Cerenkov 
cone for relativistic particles and the time distribution of the light, which is broader for 
non-relativistic particles |J , gave excellent separation between electrons and particles below 
Cerenkov threshold. 

Cosmic rays were attenuated by roughly 2 kg/ cm 2 of overburden. The cosmic ray trigger 
rate was then reduced from around 10 kHz to an acceptable level of roughly 50 Hz by an 
active veto sheild. The veto shield enclosed the detector on all sides except the bottom, as 
shown in Fig. [j] by the heavy black line surrounding the detector. The main veto shield 



14| consisted of a 15-cm layer of liquid scintillator in an external tank and 15 cm of lead 



shot in an internal tank. Following the 1993 running, additional counters were placed over 
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the crack between the endcap veto and the barrel region veto system, and below the veto 
shield along the sides. That reduced cosmic-ray background entering through veto sheild 
gaps and the bottom support structure. The combination of active and passive shielding 
tagged cosmic-ray muons that stopped in the lead shot. A veto inefficiency < 10~ 5 was 
achieved for incident charged particles. 

F. Data Acquisition 

Digitized time and pulse height of each of the 1220 detector PMTs (and each of the 292 
veto shield PMTs) were recorded when the deposited energy in the tank exceeded a threshold 
of 150 hit PMTs (~ 4MeV electron-equivalent energy) with < 4 veto PMT hits and with 
no event with > 5 veto PMT hits within the previous 15.2 /is. Activity in the detector 
or veto shield during the 51.2 jis preceding a primary trigger was also recorded, provided 
there were > 17 detector PMT hits or > 5 veto PMT hits. Data were recorded for 1 ms 
after the primary trigger at a reduced threshold of 21 PMT hits (about 0.7 MeV) in order to 
detect the 2.2 MeV 7 from neutron capture on free protons, which has a capture time of 186 
/is. The detector events were recorded without reference to the beam spill, but the state of 
the beam was recorded with the event. Approximately 94% of the recorded events occured 
between beam spills, which allowed an accurate measurement and subtraction of cosmic-ray 
background surviving the event selection criteria. 

As most muons from muon-neutrino induced events do not satisfy the PMT trigger thresh- 
old, these muons were typically past events, while the electrons from their decay were the 
primary events. In contrast, electrons from electron- neutrino induced events were usu- 
ally primary events. Future events include neutron capture 7s and (3-decay electrons and 
positrons. Identification of neutrons was accomplished through the detection of the 2.2 MeV 
7 from neutron capture on a free proton. Nitrogen and boron ground-state /5-decays oc- 
curred after the primary events with longer lifetimes of 16 and 30 ms, respectively. A given 
primary event can have many associated past events and future events. 



8 



III. NEUTRINO INTERACTIONS AND EVENT SIMULATION 



The neutrino interactions that occured in LSND came from interactions on carbon, free 
protons, and electrons in the detector liquid. All four possible neutrinos, z/ e , v e , u^, and 
Pp, contribute to neutral-current processes over the entire energy range. Charged-current 
cross sections are significantly affected by nuclear threshold effects. In the case of z/ M and 
v^ charged-current interactions, a large amount of the initial neutrino energy goes into the 
mass of the final state muon. 

Neutrino processes that are observed in LSND are classified into three categories: stan- 
dard model leptonic processes (e.g. ve — > ve elastic scattering), inverse /3-decay processes, 
and semi-leptonic processes that leave excited or fragmented nuclei in the final state. Cross 
sections in the first category may be calculated to high accuracy, better than 1%, provided 
that the neutrino energy is known. Cross sections for the inverse /3-decay reactions are in- 
ferred from the measured /?-decay lifetimes and are accurate to the order of a few percent. 
(The momentum transfers are sufficiently small that form factor dependences are well char- 
acterized.) The cross sections for the reactions involving nuclear excited states are much 
less certain ||15|| . Models such as the continuum random phase approximation (CRPA) [16j 
often require large corrections in order to account for ground state wave functions that are 
too simplistic. Fermi Gas models do not reliably take into account nuclear effects but can 
be made to produce reasonable agreement in the quasi-elastic energy region when effective 
masses are employed [118, [19 . 



We use the results of a shell model calculation ]T3| for the 12 C(v e , e~) 12 N DAR processes 
[E v < 52.8 MeV). The shell model calculation gives a similar energy shape but a lower cross 
section than the CRPA calculation |16|]. A relativistic Fermi Gas model with an effective 
mass correction employed to account for nuclear effects is used for the more energetic DIF 
neutrino processes. 

Two-body neutrino interactions are known accurately from either measurement or theory. 



Those processes are listed in Table [TV] with their associated cross section uncertainty. They 
provide the main constraints on neutrino fluxes, trigger and selection efficiencies, and other 
neutrino cross sections. Table |V| also lists the neutrino flux sources constrained by each 
of these processes. For example, the l2 C(v e ,e~) 12 A 9 . S . and the ve — > ve elastic reactions 
primarily constrain the rate of /i + DAR in the target area. Of all the 12 N states, only 
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the ground state (3 decays, and the 12 C(z/ /J , /z - ) 12 A 9 . S . reaction is the best measure of the 
7r + — > DIF rate in the target area. Those reactions that contain a final state 12 N g . s . 

have nuclear matrix elements directly related to well known nuclear matrix elements, so 
that the cross sections may be calculated to an accuracy of a few percent. The ue — > ve 
elastic reactions are Standard Model electroweak calculations and are known to better than 
a percent from the measured weak mixing angle, sin 2 9w, and the Fermi constant, Gf- 

Events in the LSND detector were simulated by using the GEANT3.21 code |17| , which 
was modified to track optical photons in addition to ionizing particles. Neutrons were 
tracked and captured on free protons via the standard MICAP interface to GEANT. Optical 
photon yields from both scintillation and Cerenkov processes were generated and tracked to 
the photomultiplier tubes. A simulation of the photomultiplier response, analog and digital 
electronics, and event trigger produced event data packets, which were a good representation 
of the LSND detector response to neutrino events. A large sample of Michel decays from 
cosmic rays was used to check the quality of simulated events. The resulting neutrino event 
samples are then processed by the same reconstruction and particle identification software 
as the beam data. The reconstructed data are then compared to beam-excess data for the 
following analyses. 



IV. DATA PROCESSING AND EVENT RECONSTRUCTION 

A new event reconstruction that improved the position resolution and the spatial cor- 
relation between the e + and neutron-capture 7 in the reaction i/ e p — > e + n was applied to 
the entire 1993-1998 data sample. Different event samples were made during the new data 
reduction, and we focus here on the measurement of electron events, which are relevant to 
the oscillation search. 

The electron selection was applied to the ~4 Terabytes of raw LSND DLT data tapes, 
using a minimal set of cuts. This process achieved roughly a 40:1 reduction in data size, while 
maintaining an 87 ± 2% efficiency for electron events, independent of electron energy above 
20 MeV. Events in this new data stream that appeared in samples from previous LSND 
analyses were labelled, and a cross-check for consistency between new and old samples was 
performed. 

Table [V| shows the electron reduction criteria and the corresponding efficiencies. First, 
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the visible energy was required to be greater than 15 MeV in order to eliminate 12 B f3 decays 
from cosmic ray /i~ that stop and capture in the oil. Second, the number of associated veto 
hits was required to be less than 4. Third, events with a laser calibration tag were rejected. 
Fourth, loose electron particle identification criteria were imposed. Fifth, the resulting 
data were subjected to a loose fiducial volume cut, which required that the reconstructed 
electron vertex be inside a volume that was greater than 10 cm from the PMT surfaces. 
Finally, cosmic-ray muon events that produced decay electrons (Michel electrons) as the 
primary event were removed. In a clean sample of cosmic-ray Michel electron events there is 
a correlation between the total number of photoelectrons at the muon time and the online 
reconstructed distance to the subsequent Michel electron; as the cosmic muon becomes more 
energetic, the distance to the Michel electron grows linearly. A two dimensional region, or 
graphical cut, was imposed to remove these events. 

The efficiency for electrons surviving the cuts was determined as follows. In an unbiased 
sample of laser-induced events with their associated accidental activities, a Monte Carlo 
(MC) electron event was inserted in place of the laser "primary". This left a MC electron 
event in the midst of the accidental events from the real laser event. The MC electrons were 
generated flat in energy and uniformly throughout the tank. Desired accidental properties of 
the laser event, e.g. veto hit count and time to activities, were preserved when the electron 
MC event was inserted. Electron reduction criteria were then applied and the efficiencies 
calculated. The removal of accidental hits in time with the laser did not significantly affect 
the efficiency measurement. 

Data were reprocessed with the new event reconstruction in order to improve the position 
resolution. The previous event reconstruction was limited due to the charge response of the 
8" PMTs used in LSND (Hamamatsu R1408). For these PMTs, the single photoelectron 
output charge distribution is approximately a broad Gaussian plus an exponential tail that 
extends well above the mean of the Gaussian. As the position and angle fits weighted the 
hit PMTs by their charge, this charge tail has the effect of smearing the reconstructed 
event positions and angles. To ameliorate this effect, a new reconstruction algorithm was 
developed that weighted the hit PMTs by a ratio of the predicted charge to the square of 
the time resolution for that predicted charge, Q pre di '^l-, and not by their measured charge, 
Qtube- (The new reconstruction also has other improvements, such as the inclusion of timing 
information in the 7 reconstruction.) This has resulted in an improvement in the position 
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correlation between the muon and the electron from stopped muon decay and between the 
neutron and the 7 from neutron capture. The mean reconstructed distance between the 
muon and decay electron improved from 22 cm with the previous reconstruction to 14 cm 
with the current reconstruction. For 2.2 MeV 7 from neutron capture, the most likely 
distance was reduced from 74 cm to 55 cm. As the accidental 7 rate is proportional to the 
cube of this distance, the resulting 7 reconstruction allows a cut on the _R 7 parameter, as 
described later in section VII, that yields a factor of two better efficiency with a factor of 
two reduction in the rate of accidental 7s. 

V. PRIMARY ELECTRON SELECTION 

The primary electron selection is next applied to the reduced data. The goal of the 
selection is to reduce the cosmic-ray background to as low a level as possible, while retaining a 
high efficiency for neutrino-induced electron events. The selection criteria and corresponding 
efficiencies are shown in Table [V|. The energy range 20 < E < 200 MeV is chosen so as to 
accept both DAR — > v e and DIF — > v e oscillation candidates. We require 20 < E e < 60 
MeV for the — > v e oscillation search and 60 < E e < 200 MeV for the z/ M — > v e oscillation 
search. Below 20 MeV there are large backgrounds from the f3 decay of 12 B created by the 
capture of stopped cosmic-ray fi~ on 12 C. Above 200 MeV the beam-related backgrounds 
from 7r + — > e + v e are large compared to any likely oscillation signal. Events with a previous 
activity within 12 /is, a future activity within 8 fis, or a bottom veto counter hit are rejected 
in order to eliminate cosmic-ray muon events. To further minimize cosmic-ray background, 
a tight electron particle identification is applied, —1.5 < x'tot < 0-5, where the allowed range 
is chosen by maximizing the selection efficiency divided by the square root of the beam- 
off background with a correlated neutron. The x'tot parameter depends on the product of 
the x parameters defined in ||. Briefly, Xr and Xa are the quantities minimized for the 
determination of the event position and direction, and Xt is the fraction of PMT hits that 
occur more than 12 ns after the fitted event time. The dependence of the x parameters 
on energy and position for Michel electrons was studied, and a correction was developed 
that made x'tot independent of energy or position. For the 1993 data only, which had an 
electronics timing problem resulting in a broader x'tot distribution, 0.3 < xlit < 0-65 was 
also required, where xfot * s computed like x'toti but with \ parameters defined in reference 
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IJ. Additionally, the trigger time is required to occur between 85 ns and 210 ns in the 500 
ns trigger window in order to reject multiple events, no veto hit is allowed within 30 ns of 
the trigger time, and the reconstructed electron vertex is required to be inside a volume 
35 cm from the faces of the photomultiplier tubes. Finally, the number of associated 7s 
with i? 7 > 10 (Ry is discussed in section VII) is required to be < 2 (< 1) for events < 60 
(> 60) MeV in order to reject neutron-induced events, which tend to have many associated 
7s. Neutrons from u e p — > e + n scattering are too low in energy (< 5 MeV) to knock out 
other neutrons; however, higher energy neutrons (> 20 MeV) typically knock out 1 or more 
neutrons. The event selection is identical for the DAR and DIF samples except for the 
associated 7 criteria. Note that the event selection is optimized for electrons in the DAR 
energy range; however, it was applied to the DIF energy range for simplicity and in order 
that a common selection criteria be used over the entire interval from 20 — 200 MeV for 
oscillations from both DAR u^ — > v e and DIF u^ — > u e . 

In addition to the electron reduction and selection efficiencies, Table |V] also shows the 
efficiencies due to the data acquisition (DAQ) and veto deadtime. The total efficiency for 
electrons in the fiducial volume with energies in the range 20 < E e < 60 MeV is 0.42 ±0.03. 

VI. CONVENTIONAL NEUTRINO PROCESSES 

The neutrino oscillation analysis consists of two steps. The first step is to determine 
the best values for the numbers of events from standard neutrino processes in a way that 
minimizes the systematic uncertainty due to the electron selection. The second step is to use 
those measured neutrino backgrounds as central values in a fit to the oscillation parameters, 
allowing the backgrounds to vary around the central values within their uncertainty. The 
first step will be discussed in this section. 

The inclusive electron data set provides a common selection for all neutrino processes 
important to the oscillation analysis. Some of these have well-determined cross sections: the 
12 N ground state events, ue elastic events, and up charged- current events. These events serve 
to constrain the neutrino fluxes and the selection efficiencies. They also provide important 
constraints on uncertain cross sections, such as 12 N excited state events, where the nuclear 
response function is not well known. 

Once the primary electron is selected, events are categorized by whether or not there are 
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associated events in the past or future of the primary. This categorization isolates most of 
the important reactions. The simplest event topology has an electron with no correlated 
event in the past or future, i.e. inclusive electrons. Ground state events are selected by 
searching for 12 N (3 decay within 70 ms and 70 cm of the primary electron event. Events 
from the process u e p — > e + n have a correlated 7 from neutron capture within 1 ms. Muon 
neutrino induced events are selected efficiently because the Michel electron decay of the 
muon satisfies the primary electron criteria. The additional requirement of a prior event 
within 10 /is efficiently finds the initial neutrino interaction muon event. The muon events 
can have, in addition, future events from neutron-capture 7s and nuclear j3 decays. A list 



of the various event categories is shown in Table VI 



A least squares fit was designed to find the best values for the neutrino fluxes, efficiencies, 
and cross sections. It fits those parameters by minimizing the \ 2 formed from the predicted 
number of events in various distributions compared to the observed number of events. The 



distributions are chosen to be sensitive to each of the parameters in question. Table [VTl] lists 
the parameters adjusted in the fit, along with the fitted correction values, central correction 
values, and nominal parameter values. The central correction value is the Gaussian error 
by which each parameter was allowed to vary in the fit. The final fitted value for each 
parameter is the product of the nominal parameter value and the fitted correction value. 
The agreement between the data and the least squares fit is good. The fitted DIF neutrino 
flux and tt~/tt + ratio are about one sigma lower than the nominal values; however, the 
nominal values are used when estimating the neutrino background to the oscillation search. 
The v^C — > fi~N* cross section is lower than current theoretical predictions |15|], [lf| but is 
in agreement with our earlier measurement [^j] . 

Fig. [5] shows the electron and j3 energy distributions and the time between the electron 
and {3, At, for 12 C(z/ e , e~) 12 N g s . scattering events. The energy and angular distributions 
for inclusive electron events are shown in Fig. |^, where E e is the electron energy and 9 V is 
the angle between the incident neutrino and outgoing electron directions. Neutrino-electron 
elastic scattering events are clearly visible near cos6 u ~ 1. Fig. [7| shows the angular 
distribution in more detail (top plot) as well as the energy distribution (bottom plot) for 
the neutrino-electron elastic scattering events with cos8 u > 0.9 and with 12 C{v ei e~) 12 N gs ^ 
events removed. Fig. |8] shows the muon and (3 energy distributions and the time between 
the muon and /3 for 12 C(^, fi~) 12 N g . s . scattering events. Finally, Fig. |9| shows the Michel 



14 



electron and muon energy distributions, the time between the muon and electron, At, and 
the distance between the reconstructed electron and muon positions, Ar, for v^C — > n~N, 
v^C — > fi + B, and v^p — > fi + n inclusive scattering events. Cross sections for uC scattering 
based on a partial data sample have been published previously P0|, [21]. Final cross sections 
for ve elastic scattering v e C scattering |2"3(] , and u^C scattering [|24J will be reported 
elsewhere and are consistent with the nominal parameter values shown in Table |V11| . 



VII. THE DECAY-AT-REST OSCILLATION ANALYSIS 



A. Signal and Background Reactions 

The primary oscillation search in LSND is for — > v e oscillations, where the arise from 
/i + DAR in the beam stop and the v e are identified through the reaction v e p — > e + n. This 
reaction allows a two-fold signature of a positron with a 52 MeV endpoint and a correlated 
2.2 MeV 7 from neutron capture on a free proton. There are only two significant neutrino 
backgrounds with a positron/electron and a correlated neutron. The first background is from 
/i~ DAR in the beam stop followed by v e p — > e + n scattering in the detector. As mentioned 
earlier, this background is highly suppressed due to the requirements that a n~ be produced, 
the 7r~ decays in flight, and the \x~ decays at rest prior to capture. The v e flux is calculated 
to be only ~ 8 x 10~ 4 relative to the flux in the 20 < E v < 52.8 MeV energy range. The 
second background is from 7r~ DIF in the beam stop followed by u^p — > fi + n scattering in the 
detector. (Additional contributions are from v^C — > fi + nX and v^C — > fi~nX scattering.) 
This background will mimic the oscillation reaction if the fi + is sufficiently low in energy 
that it is below the threshold of 18 hit PMTs, corresponding to < 4 MeV. Table jVTJ I 
shows the estimated number of events in the 20 < E e < 60 MeV energy range satisfying 
the electron selection criteria for 100% v e transmutation and for the two beam-related 
backgrounds with neutrons. Uncertainties in the efficiency, cross section, and v flux lead to 
systematic errors of between 10% and 50% for the signal and backgrounds discussed below. 

The largest beam-related background with a correlated neutron is due to v e produced 
in the beam stop by conventional processes. Such events are identical to the oscillation 
candidates, and are identified via the reaction v e p — > e + n. Their most important source is 
the DAR of /i~ in the beam stop. The total background due to intrinsic v e in the beam is 
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the product of neutrino flux (1.08 x 10 n z/ e /cm 2 ), average cross section over the entire energy 
range (0.72 x 10~ 40 cm 2 ) ||25|| , the number of free protons in the fiducial volume (7.4 x 10 30 ), 
the fraction of events with E > 20 MeV (0.806), and the average positron reconstruction 
efficiency after cuts (0.42), which gives a total background of 19.5 ± 3.9 events before any 7 
selection. Another possible source of z/ e , the direct decay of 7r~ — > e~u e , is negligible, as a 
consequence of its low branching ratio (1.2 x 10 -4 ), the 1/8 ratio of tt~ to 7r + in the target, 
and the capture of tt~ in the material of the beam dump. 

A related background is due to v e 12 C — > e + 11 B n scattering. The cross section to 
the 12 B ground state is calculated to be 6.3 x 10~ 42 cm 2 p6| , and the cross section to the 
11 B n final state is estimated to be at least a factor of two smaller, especially because the 
first four excited states of 12 B are stable against neutron emission. Therefore, we estimate 
that this background is < 2% of the v e p —>■ e + n background and is negligible. Furthermore, 
the maximum positron energy from this background is 36.1 MeV, so that almost all of the 
positrons are below 36 MeV. 

The second most important source of beam-related background events with correlated 
neutrons is the misidentification of and charged- current interactions as v e events. 
Because of the energy needed to produce a fi, such a or must arise from a 7r that decays 
in flight. In the tank the interacts by either v^p — > fi + n or (less often) z/^C — > [i + nX, 
followed by fi + — > e + v e i>^. The interacts by u^C — > fi~nX. There are four possible 
reasons for the misidentification. First, the muon can be missed because the fi + lifetime is 
> 12/xs or the deposited energy is below the 18 phototube threshold for activity triggers. 
The latter can occur either because the muon is too low in energy or is produced behind the 
phototube surfaces. The detector Monte Carlo simulation is used to show that this threshold 
corresponds to a fi kinetic energy, T M , of approximately 3 MeV. The background rate from 
u^p — > fi + n is written as the product of the total flux above threshold (2.56 x 10 11 ^/cm 2 ), 
the average flux-weighted cross section (4.9 x 10~ 40 cm 2 ) J25j, the fraction of /i + having < 
3 MeV or r M > 12/is (0.0258), the number of free protons in the fiducial volume (7.4 x 10 30 ), 
the positron efficiency (0.42), and the fraction of events with E > 20 MeV (0.816), for a 
background of 8.2 events. Similar estimates for the backgrounds from v^C — > fi + nX and 
u^G — > iTnX |27j add 0.4 and 1.4 events, respectively, for a total of 10.0 ± 4.6 events. It 
is estimated [27] that about 80% of the v^C — > fi + X and 6% of the v^C — > fi~X scattering 
events will have a recoil neutron. 



16 



Second, a fi above the hit threshold can be missed if a prompt decay to e caused the muon 
and electron to be collected in a single event which is then misidentified as an e. This effect is 
considerably suppressed by the electron selection and the requirement that the reconstructed 
time be consistent with the triggered event time. The detector Monte Carlo simulation 
shows that this misidentification only occurs for /i + decays within 100 ns, decreases with T M 
, and is almost zero above 10 MeV. Using the Monte Carlo misidentification probabilities, a 
calculation similar to that above implies a background of 0.2 ±0.1 events. 

Third, the /i + can be lost because it is produced behind the PMT surface and the electron 
radiates a hard 7 that reconstructs within the fiducial volume. A background of 0.2 ± 0.1 
events is estimated from the Monte Carlo simulation. 

Fourth, a muon can be missed by trigger inefficiency After 1994, we acquired for many 
online positron triggers complete digitization information for all veto and detector PMTs 
over the 6 [is interval prior to the positron. Analysis of these data, discussed below, shows 
the trigger inefficiency for low-energy muons to be negligible. 

There are additional backgrounds from v e produced by [T — > e~z/ M z/ e and 7r~ — > e~v e 
DIF. These v e can interact on either Cora free proton to yield the oscillation signature of 
a positron and a recoil neutron. For 20 < E e < 60 MeV, 0.1 ±0.1 events are estimated. The 
reactions u e 12 C — > e~~nX and v e 13 C — * e~nX are negligible (< 0.1 events) over the 20 < 
E e < 60 energy range and cannot occur for E e > 20 MeV and E e > 36 MeV, respectively. 
Other backgrounds, for example v^C — ► u^nyX with E n > 20 MeV, v e C — > e~pX followed 
by 13 C(p, n) 13 N, and v^C — > n~X followed by \T~ capture, are also negligible. 

The total background due to pion and muon DIF is 10.5±4.6 events before any 7 selection. 
It has a detected energy spectrum which is very close to that for positrons from /i + decay. 

A final source of background is neutrons from the target that find their way into the 
detector tank. However, a stringent limit on beam neutron background relative to the 
cosmic neutron background has been set by looking for a beam-on minus beam-off excess of 
neutron events that pass neutron PID criteria in the 40-180 MeV electron equivalent range 
0. No excess has been observed, which implies that the beam-related neutron background 
is less than 1% of the total beam-unrelated background and is negligible. 

The number of events expected for 100% — > i> e transmutation followed by u e p — > 
e + n scattering (plus a small contribution from v e C — > e + Bn scattering) is 33300 ± 3300 
events, where the systematic error arises from uncertainties in the neutrino flux (7%) and 
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e + efficiency (7%). This number of events is the product of neutrino flux (1.26 x 10 14 z//cm 2 ), 
the average cross section over the entire energy range (0.95 x 10~ 40 cm 2 ), the average 
positron reconstruction efficiency (0.42), the fraction of events with E > 20 MeV (0.894), 
and the number of free protons in the fiducial volume (7.4 x 10 30 ). 



B. The Positron Criteria 



The positron/electron selection criteria (LSND is insensitive to the sign of the charge) 
for this primary oscillation search is described in detail in section V. 



C. The Correlated 2.2 MeV 7 Criteria 

Correlated 2.2 MeV 7 from neutron capture are distinguished from accidental 7 from 
radioactivity by use of the likelihood ratio, i? 7 , which is defined to be the likelihood that 
the 7 is correlated divided by the likelihood that the 7 is accidental. i? 7 depends on three 
quantities: the number of hit PMTs associated with the 7 (the multiplicity is proportional 
to the 7 energy), the distance between the reconstructed 7 position and positron position, 
and the time interval between the 7 and positron (neutrons have a capture time in mineral 



oil of 186 /is, while the accidental 7 are uniform in time). Fig. 10 shows these distribu- 
tions, which are obtained from fits to the data, for both correlated 2.2 MeV 7 (solid curves) 
and accidental 7 (dashed curves). To determine i? 7 , the product of probabilities for the 
correlated distributions is formed and divided by the product of probabilities for the un- 
correlated distributions. The accidental 7 efficiencies are measured from the laser-induced 
calibration events, while the correlated 7 efficiencies are determined from the Monte Carlo 
simulation of the experiment. Similar results for the correlated 7 efficiencies are obtained 
from the cosmic-ray neutron events, whose high energy gives them a slightly broader po- 



sition distribution. The efficiencies for different i? 7 selections are shown in Table [LX|. The 



systematic uncertainty of these efficiencies is estimated to be ±7% of their values. Note that 
with the new reconstruction, the correlated 7 efficiency has increased while the accidental 7 
efficiency has decreased. For i? 7 > 10, the correlated and accidental efficiencies are 0.39 and 
0.003, respectively. For the previous reconstruction the R° > 30 cut gave correlated 
and accidental efficiencies of 0.23 and 0.006, respectively. 
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As checks of the likelihood distributions, Fig. |TT| shows the _R 7 distributions for v e C — > 
e~~Ng. s . exclusive events [p0[ , where the N g , s . (3 decays. By definition, the v e C — > e~N g , s , 
reaction has no recoil neutron, so that its Ry distribution should be consistent with a purely 
accidental 7 distribution. A fit to the Ry distribution finds that the fraction of events with a 
correlated 7, f c , is f c = -0.004 ±0.007 (x 2 = 4.6/9 DOF). Fig. shows the Ry distribution 
for the sample of /i events arising from the reactions v^C — > fi~X, v^C — > fi + X, and 
v^p — * fi + n. Correlated 7 are expected for ~ 14% of these events [p7| . A fit to the Ry 
distribution gives f c = 0.129 ± 0.013 (x 2 = 8.2/9 DOF), in agreement with expectations. 



Fig. [13] shows the distributions of Ar, At, and Nhits for events with Ry > 1 (left side) 
and R 7 < 1 (right side). The top plots show the distance between the reconstructed 7 
position and positron position, Ar, the middle plots show the time interval between the 7 
and positron, At, and the bottom plots show the number of hit PMTs associated with the 
1, N hUs . 



D. Neutrino Oscillation Results 

Table |X| shows the statistics for events that satisfy the selection criteria for the primary 
V/j, — > v e oscillation search. An excess of events is observed over that expected from beam-off 
and neutrino background that is consistent with neutrino oscillations. A x 2 fit to the i? 7 
distribution, as shown in Fig. |l|, gives f c = 0.0567 ± 0.0108 (x 2 = 10.7/9 DOF), which 
leads to a beam on-off excess of 117.9 ± 22.4 events with a correlated neutron. Subtracting 
the neutrino background from /i~ DAR followed by u e p — > e + n scattering (19.5 ±3.9 events) 
and 7r~ DIF followed by u^p — * fi + n scattering (10.5 ±4.6 events) |28[ leads to a total excess 
of 87.9 ± 22.4 ± 6.0 events, as shown in Table [XI|. This excess corresponds to an oscillation 
probability of (0.264 ± 0.067 ± 0.045)%, where the first error is statistical and the second 
error is the systematic error arising from uncertainties in the backgrounds, neutrino flux 
(7%), e + efficiency (7%), and 7 efficiency (7%). Note that our previously published result 
§, based on the 1993-1995 data sample, was (0.31 ± 0.12 ± 0.05)%. Table [XTj shows the 
effect on the fitted oscillation probability of tightening some of the selection criteria. 

A clean sample of oscillation candidate events can be obtained by requiring Ry > 10, 
where as shown in Table |X|, the beam on-off excess is 49.1 ± 9.4 events while the estimated 



neutrino background is only 16.9 ± 2.3 events. Fig. 15 shows the individual 7 distributions 



19 



for events with 20 < E e < 60 MeV and with i? 7 > 1 (left side) and R y < 1 (right side). The 
top plots show the distance between the reconstructed 7 position and positron position, Ar, 
the middle plots show the time interval between the 7 and positron, At, and the bottom 
plots show the number of hit PMTs associated with the 7, Nh its . Fig. ^ displays the energy 
distribution of events with i? 7 > 10. The shaded regions show the combination of neutrino 
background plus neutrino oscillations at low Am 2 . The data agree well with the oscillation 
hypothesis. As mentioned in section I, the 1993-1995 data runs employed a 30 cm water 
target, while the 1996-1998 data runs used a high-Z metal target. A comparison of the 



energy distributions of the two data samples is displayed in Fig. [17], which shows that the 
data samples are consistent within statistics. 

Fig. [18] shows the cos Q v distribution for events with i? 7 > 1 and 36 < E e < 60 MeV. This 
energy range is chosen because it is particularly clean with reduced neutrino background, so 
that the u e p —> e + n reaction should dominate, while the 7 requirement is relaxed to increase 
the statistics. 9 is the angle between the incident neutrino and outgoing positron directions. 



The shaded region in Fig. [18] shows the expected distribution from a combination of neutrino 
background plus neutrino oscillations. The < cosO^ >= 0.04 ± 0.12, in agreement with the 
expectation of ~ 0.12. 

Figs. [Tj^ (D > 10 cm ) and (D > 35 cm ) show the spatial distributions for events with 
R y > 10 and 20 < E e < 60 MeV, where z is along the axis of the tank (and approximately 
along the beam direction), y is vertical, and x is transverse. The shaded regions in Figs. 



19| and |20] show the expected distributions from a combination of neutrino background plus 



neutrino oscillations. Fig. ^T] shows scatter plots of the x-y and y-z spatial distributions 



for events with _R 7 > 10, 20 < E e < 60 MeV, and D > 35 cm. Figs. || and |3] show 
the x'tot an d veto hit distributions for events with i? 7 > 10 and 20 < E e < 60 MeV. The 
solid histogram in the veto hit figure shows the distribution from v e C —> e~N g s scattering. 



Finally, Fig. |24] shows the L v j 'E v distribution for events with i? 7 > 10 and 20 < E e < 60 
MeV, where L v is the distance travelled by the neutrino in meters and E u is the neutrino 
energy in MeV determined from the measured positron energy and angle with respect to the 
neutrino beam. The data agree well with the expectation from neutrino background plus 
neutrino oscillations at low Am 2 (x 2 = 4.9/8 D.O.F.) or high Am 2 (x 2 = 5.8/8 D.O.F.). 
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E. Tests of the — > v e Oscillation Hypothesis 

A variety of tests of the — > v e oscillation hypothesis have been performed. One test 
of the oscillation hypothesis is to check whether there is an excess of events with more than 
one correlated 7. If the excess of events is indeed due to the reaction v e p — > e + n, then 
there should be no excess with more than one correlated 7 because the recoil n is too low 
in energy (< 5 MeV) to knock out additional neutrons. If, on the other hand, the excess 
involves higher energy neutrons (> 20 MeV) from cosmic rays or the beam, then one would 
expect a large excess with > 1 correlated 7, as observed in the beam-off cosmic ray data. 
However, as shown in Table |X11J|, the excess of events with more than one correlated 7 



is approximately zero for both the full 20 < E e < 60 MeV energy region and the lower 
background 36 < E e < 60 MeV energy region, as expected for the reaction v e p — > e + n. 

Another test of the oscillation hypothesis is to check the "event lookback" for events 
that satisfy the oscillation criteria in order to ensure that the v^p — > fi + n background is 
calculated correctly. The "event lookback" was installed prior to the 1995 running and 
consisted of an extra trigger that read out all hit detector PMTs in the 6 fis interval before 
a primary event. Any background just below the 18 hit muon threshold will clearly show 
up, especially in the hit range with > 11 lookback hits, where the probability of having an 
accidental lookback is only 5.6%. However, as shown in Table [XI V| , the excess of primary 
events with i? 7 > or i? 7 > 10 is consistent with the accidental lookback probability. Thus, 
this "event lookback" check provides additional assurance that the v^p —> fi + n background 
calculation of 10.5 events is not underestimated. 

Both major backgrounds with a correlated neutron arise initially from 7r~ DIF. Therefore, 
a final test of the oscillation hypothesis is to check whether the flux from tc~ DIF is correct. 
However, this has already been tested by the fit to the i? 7 distribution, discussed above, 
of VyC — > /Ji~X, v^C — > /i + A, and v^p — > fi + n inclusive events ||21|| . For these reactions, 
correlated 7 are expected for ~ 14% of the events p7| , due mainly to the flux. A fit to the 
R y distribution gives f c = 0.129 ± 0.013 (y 2 = 8.2/9 DOF), in agreement with expectation. 
Fig. p5| shows the time to the previous event for i? 7 > 10 electron events prior to applying the 
At pas t > 12/xs selection. The top graph in the figure shows that the beam excess events are 
in agreement with our expectations for the v^p — > fi + n and 12 C(z/ /J , 12 B* channels. With 
the same data on a smaller vertical scale, the bottom graph shows events with accidental 
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past activities, in agreement with expectations from random cosmic ray backgrounds. 



VIII. THE DECAY-IN-FLIGHT OSCILLATION ANALYSIS 

The high energy v^s from tt + decay-in-flight are a potential source of — > v e oscillation 
events. The contamination of z/ e s from standard sources is small, at the level of 0.1%. 
Unfortunately, the cosmic-ray backgrounds are large, with the result that the fluctuations 
in the beam-on, beam-off subtraction are comparable to the expected signal. Prior to 1996, 
it was realized that the 1996-1998 data would not support a stand alone, decay-in-flight 
analysis due to the larger beam-off backgrounds that are inherent in running with a heavy 
target. However, the analysis presented here is extended up to an electron energy of 200 
MeV because the decay-in-flight data constrain the region > 2 eV 2 , especially around 6eV 2 . 

The above analysis is applied to data in the energy range 60 < E e < 200 MeV , with 
the additional requirement that there be no associated 7. This sample is only sensitive 
to v^ — > v e oscillations, and results in a beam on-off excess of 14.7 ± 12.2 events. The 
signal expected for 100% v^ to v e transmutation is estimated to be 7800 events, and the 
v e background from fi + — > e + v fl v e , 7r + —> e + v e , and ve — > ve is estimated to be 6.6 ± 1.7 
events, resulting in a total excess of 8.1 ± 12.2 ± 1.7 events or an oscillation probability of 



(0.10 ± 0.16 ± 0.04)%, as shown in Table [XV| . This result is lower than but consistent with 



our higher precision analysis of the 1993-1995 data sample 0. That analysis determined 
the selection criteria by maximizing the acceptance divided by the square root of the beam- 
off background, which produced much less beam-off background overall. It gave a total 
excess of 18.1 ± 6.6 ± 4.0 oscillation events, corresponding to an oscillation probability of 
(0.26±0.10±0.05)%. Due to changes in the neutrino production target, the 1996-1998 data 
sample had reduced DIF flux and higher beam-off background compared to the 1993-1995 
data. Based on our DAR oscillation result and assuming that CP is conserved in the lepton 
sector, we would expect the DIF oscillation probability to be ~ 0.26% at high Am 2 (where 
(1.27 Am 2 L u /E v ) » 1) and ~ 0.05% at low Am 2 (where (1.27Am 2 L u /E u ) << 1). 

The 0.16% statistical error on the oscillation probability in the present analysis is larger 
than the 0.10% statistical error of the previous analysis. That is because the present analysis 
uses the electron selection criterion developed for the DAR region below 60 MeV, which is 
less effective in removing the background to electron events in the DIF region above 60 
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MeV. When compared to the previously published DIF analysis, the beam-off background 
for these higher energy events is 3.42 times higher in the present analysis, while the number 
of expected events for 100% — > v e transmutation is only 1.16 times higher. The previous 
analysis observed a 2.6 sigma excess, compared to the 0.6 sigma excess of the present analysis. 

IX. THE NEUTRINO OSCILLATION FIT 

A. Introduction 

We describe in this section the (sin 2 29, Am 2 ) likelihood (£) fitter. The fitter is applied to 
beam-on events in the final oscillation sample and calculates a likelihood in the (sin 2 29, Am 2 ) 
plane in order to extract the favored oscillation parameters. The fit is similar to that 
performed in reference |[29|| . 

The £ product in the (sin 2 29, Am 2 ) plane is formed over the individual beam-on events 
that pass the oscillation cuts. This three-dimensional contour is sliced to arrive finally at the 
LSND allowed oscillation region. The beam-related backgrounds are determined from MC 
event samples for each individual background contribution. The MC contains the trigger 
simulation and generally very well reproduces the tank response to all particles of interest. 
Agreement between the data and MC is excellent. The fit is over the entire electron energy 
range 20 < E e < 200 MeV. Therefore, DIF oscillations and DIF backgrounds in addition to 
the usual DAR processes are considered. 

B. Formalism 

Each beam-on event is characterized by four variables: the electron energy, E e , the 
electron reconstructed distance along the tank axis, z, the reconstructed direction cosine 
the electron makes with the neutrino, cos^, and the likelihood ratio that the event has a 
correlated 2.2 MeV 7, i? 7 . Each of the neutrino-induced background processes is simulated, 
and the simulation is compared to real events in the detector. Accidental 7 events are used 
with real neutrino processes to simulate accidental events. Beam-off events are used as 
a background contribution after scaling by the measured time-dependent duty factor. The 
duty factor for this analysis was determined by using the entire raw event sample to measure 
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the ratio of beam-on time to beam-off time. The raw event sample consists mostly of beta- 
decay events and is, to a good approximation, unbiased by beam-related events. The duty 
factor for each run was determined by dividing the number of raw events when the beam-on 
bit was set by the number of raw events where it was not set. This resulted in a duty factor 
for each run that was used to weight the beam-off events to determine the beam-unrelated 
subtraction for the final event sample. 

For every point in the (sin 2 29, Am 2 ) plane, oscillation signal events are generated to 
complete the description of sources expected in the beam-on sample. There are 5697 beam- 
on events in the data sample, and a likelihood is calculated for each one based on the values 
of E e , z, cos9 u and R 1 . 

Formally, each neutrino beam-on event j is assigned a probability Pj(E e , z, cos9 u , i? 7 ) 
equal to a sum of probabilities qi(E e , z, cos V , Ry) from the backgrounds plus oscillations. 
It then remains to add the % with expected fractional contributions r\ and take the product 
over all the beam-on events. The likelihood is thus 

N beam — on 

£=( n Pi), (9-i) 

where 

N contributions 

Pj(E ej , R 7j , cos 9 vj , Zj) = <li( E ej, Ryj, cos 9 vj , Zj) ■ r { . (9.2) 

i=i 

Additionally, two normalization requirements must hold: 

N contributions 

E r i = l, (9.3) 

i=i 

and 

J dE e dRyd(cosd u ) dz qi(E e , R 7 ,cos9 u , z) = 1 (9.4) 

for each contribution, i. Together, these requirements ensure that every observed beam-on 
event has a probability of occurrence equal to 1. 



C. Background Variation 

It is necessary to allow for the fact that the backgrounds are not perfectly known. The 
background variation is performed by calculating the above likelihood at each point in the 
(sin 2 29, Am 2 ) plane many times, varying over the expected a for each background. For each 
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background configuration, the C is weighted with a Gaussian factor for each background that 
is off its central value. The background configurations are varied so that the beam-unrelated 
background (BUB) varies independently and the beam-related backgrounds (BRBs) are 
locked together. Different background varying procedures give very similar results. 

D. The expression for the Likelihood 

Finally, the likelihood can be expressed as 



where the / VN bgd represents, schematically, the background variation described above. 

E. The Input 

The qi for each of the background and signal processes are all generated from the MC, 
except for the BUB q^, which is generated from the beam-off data events. There are separate 
MC runs for each of the above BRB processes. Some of these backgrounds are grouped 
together (appropriately weighted) into a few common gjS for easier bookkeeping, as indicated 
in Table |XV1| . This is done for backgrounds which don't need to be separately varied. 
Several small, beam-related backgrounds, DIF ve — > ue elastic scattering and 7r + — > e + u e 
DAR followed by v e C — > e~N scattering, are contained in their DAR and DIF counterparts. 

F. Slicing the contour 

1. The Feldman- Cousins Method 

The Feldman-Cousins method |30| can be applied to the LSND C contour in the following 
way. At a particular point in the (sin 2 29, Am 2 ) plane, create thousands of generated data 
sets comprised of background and oscillations. For each Monte Carlo experiment compute 
5L = L Max — L MC , where L = logC, L MC is L at the particular point in the (sin 2 29, Am 2 ) 
plane assumed in the Monte Carlo, and L^ax is the log likelihood at the values of sin 2 29 
and Am 2 that maximize C From a histogram of SL for the thousands of Monte Carlo data 
sets one obtains the selection that contains, for example, 90% of the experiments. Finally, 




i=l 



(9.5) 



25 



determine this selection at many points in the (sin 2 29, Am 2 ) plane. The resulting function 
of sin 2 29 and Am 2 corresponds to the 90% CL. allowed LSND region. 

This approach, as practiced in reference |29| , required large amounts of CPU. Even scan- 
ning a judiciously chosen (sin 2 29, Am 2 ) region is CPU intensive, and setting up and run- 
ning the generated data sets would take many months. Therefore, the full Feldman- Cousins 
method will not be followed here. As shown below, using slices derived from a different 
LSND data set to determine the L contours for this data set, the results obtained with the 
Feldman-Cousins method are similar to other methods. 



2. The Bayes Method 

For the Bayes method one presumes a prior expectation of the oscillation parameters 
from 0.01 to 100.0 eV 2 in Am 2 and 0.001 to 1.0 in sin 2 29. The assumption of this prior 
expectation is what makes this approach Bayesian. Each bin in the (sin 2 29, Am 2 ) plane is 
assigned a weight w, where w = 5x 5y-C That is, the weight is the measure of the probability 
distribution times the C The measure SxSy is taken to be <5(lnsin 2 29)6 {hi Am 2 ). The 90% 
and 99% C.L. regions are then determined by integrating over the (sin 2 29, Am 2 ) plane. 

3. The Constant- Slice Method 

The constant-slice method makes a slice at a constant value of L. If, for example, the log 
likelihood were a two-dimensional Gaussian, slices of 2.3 and 4.6 units down from the peak 
L would correspond to 90% and 99% C.L., respectively. Fig. shows that the Feldman- 
Cousins, Bayesian, and constant-slice methods all give about the same 90% regions. Note 
that for the Feldman-Cousins method the slices are derived from a different LSND data 
set. We use the constant-slice method in this paper to denote the favored regions in the 
(sin 2 29, Am 2 ) plane. 

G. Statistical Issues and Technical Hurdles 

Preserving correlations in the E e ,R 1 ,cos9 v ,z parameter space over which the C fit is 
performed is sometimes difficult, due to the fact that for certain backgrounds the 3600 bin 
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parameter space is too large to characterize. In particular, i? 7 , with its logarithmic behavior 
for backgrounds in which uncorrelated 7s are present, is especially difficult. This problem 
was resolved for the BUB by binning the other parameters very coarsely, effectively ignoring 
correlations in some regions of the four- dimensional parameter space. Statistical problems 
with the MC BRB sample, in which uncorrelated 7s are present, were dealt with in a similar 
manner. Such measures were safe approximations for the fiducial volume of interest. 

Other technical difficulties in certain ranges of Am 2 were overcome with weighting tech- 
niques. The origin of the difficulties was always one of limited statistical samples that 
characterized the probability distribution functions for the backgrounds. Another problem 
involved re- weighting by sin 2 (1.27Am 2 ^), which required prohibitive numbers of MC events 
and the simultaneous breaking of correlations in the four- dimensional space. However, these 
difficulties were overcome by smearing L u , the distance travelled by the neutrino, and E u , 
the neutrino energy, with the Gaussian widths determined from the position and energy 
resolutions. 



H. Results 



A (sin 2 29, Am 2 ) oscillation parameter fit for the entire data sample, 20 < E e < 200 



MeV, is shown in Fig. \27[ The fit includes both — > v e and — > v e oscillations, as 
well as all known neutrino backgrounds. The inner and outer regions correspond to 90% 
and 99% CL allowed regions, while the curves are 90% CL limits from the Bugey reactor 



experiment , the CCFR experiment at Fermilab |32[] , the NOMAD experiment at CERN 



3~3fl , and the KARMEN experiment at ISIS [34|. The most favored allowed region is the 



band from 0.2 - 2.0 eV 2 , although a region around 7 eV 2 is also possible, but has been made 
less probable by the — > v e analysis. 

The KARMEN experiment also searches for — > i> e oscillations with a detector that 
is similar to LSND. A comparison of the two experiments is given in Table [XV11| . LSND 
is a more massive detector, has a higher intensity neutrino source, and has good particle 
identification, while KARMEN has better energy resolution and the advantage of a much 
lower duty factor that helps eliminate cosmic-ray events. In addition, KARMEN is located 
17.5 m from the neutrino source, compared with 30 m for LSND. Therefore, the experiments 
have sensitivities that peak at different values of Am 2 . At low Am 2 , for example, an 
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experiment at 30 m is 2.94 times more sensitive to neutrino oscillations than an experiment 
at 17.5 m. Note that a global analysis of the two experiments was performed by Eitel |29| 
using intermediate data sets. 

The event breakdown from the 20 < E < 200 MeV four-dimensional fit is shown in Table 
XV1| at the best-fit point 



(sin 2 29, Am 2 ) best . fit = (0.003, 1.2eV 2 ). 



The number of z/ M — * v e oscillation events at the best-fit point is 89.5 events, which agrees 
well with the 87.9 ± 22.4 ± 6.0 event excess from the fit to the _R 7 distribution. The whole 
low Am 2 region gives an almost equally good fit within 0.5 log-likelihood units. Projections 



onto E e , R^, z, cos 9 V from the four-dimensional fit at the best fit value of (sin 2 26, Am 2 ) are 



plotted in Fig. 28. The fit is relatively insensitive to the starting values and gives good 



overall agreement wth the data. 



X. CONCLUSIONS 



The final LSND — > v e oscillation results are presented for all six years of data collection, 
1993-1998. The analysis employed a new event reconstruction that greatly improved the 
correlation of the e + and 2.2 MeV 7 from the reaction v e p — > e + n, thus greatly reducing 
the background from neutrino events followed by an accidental 7. These final results are 
consistent with our earlier analysis of the 1993-1995 data sample 0; in particular, the results 
from the 1993-1995 data sample, which used a water target, are consistent with the results 
from 1996-1998, which made use of a high-Z target. 

A global fit was performed to all event categories shown in Table [V| in order to check 
our understanding of neutrino processes in the experiment. The parameters resulting from 
this fit, shown in Table |V1J| , together with neutrino oscillations, yield a good description of 
all the observed data. 

The LSND experiment provides evidence for neutrino oscillations from the primary — > 
i> e oscillation search. A total excess of 87.9 ± 22.4 ± 6.0 v e p — > e + n events with e + energy 
between 20 and 60 MeV is observed above expected neutrino-induced backgrounds. This 
excess corresponds to an oscillation probability of (0.264 ± 0.067 ± 0.045)%. A fit to all 
of the LSND neutrino processes determines the allowed oscillation parameters in a two- 
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generation model. In conjunction with other available neutrino oscillation limits, the LSND 
data suggest that neutrino flavor oscillations occur with a Am 2 in the range 0.2 — 10 eV 2 /c 4 . 

In addition, using the same event selection, results are also presented for the decay-in- 
flight energy region. Although a clear event excess is not observed, the results are consistent 
with the Dp — *> v e oscillation signal and with our higher precision analysis of the 1993-1995 
data sample [[J, which determined the selection parameters by maximizing the acceptance 
divided by the square root of the beam-off background and which had much less beam-off 
background overall. 

At present, the LSND results remains the only evidence for appearance neutrino oscil- 
lations and implies that at least one neutrino has a mass greater than 0.4 eV/c 2 . The 



MiniBooNE experiment at Fermilab |5S] , which is presently under construction, is expected 
to provide a definitive test of the LSND results, and if the neutrino oscillation results are 
confirmed, will make a precision measurement of the oscillation parameters. 
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TABLE I: The proton beam statistics for each of the years of running, 1993 through 1998. 



Year 


Charge (C) 


Protons (xlO 22 ) 


A6 target 


Active Targets 


1993 


1787 


1.12 


water 


Al, A2, A6 


1994 


5904 


3.69 


water 


Al, A2, A6 


1995 


7081 


4.42 


water 


Al, A2, A6 


1996 


3789 


2.37 


high-Z metal 


A6 & partial A2 


1997 


7181 


4.48 


high-Z metal 


A6 only 


1998 


3154 


1.97 


high-Z metal 


A6 only 



TABLE II: The absorption rates for materials in the target area [12] 



Material 


z 


/i Absorption Rate ([is ) 


H 


1 


0.00042 ± 0.00002 


Be 


4 


0.0074 ± 0.0005 


C 


6 


0.0388 ± 0.0005 





8 


0.1026 ±0.0006 


Al 


13 


0.7054 ± 0.0013 


Fe 


26 


4.411 ±0.024 


Cu 


29 


5.676 ± 0.037 


Zn 


30 


5.834 ± 0.039 


Mo 


42 


9.61 ±0.15 


Ta 


73 


12.86 ±0.13 


Pb 


82 


13.45 ±0.18 


U 


92 


12.60 ±0.04 
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TABLE III: Average neutrino fluxes in LSND. Both decay at rest (DAR) and decay in flight (DIF) 
are shown in v/cm 2 . The v^ and v^ DIF fluxes are above [i production threshold. 



Source 


Type 


1993-1995 Flux 


1996-1998 Flux 


Total Flux 


H+ DAR 


v^ and v e 


7.38 x 10 13 


5.18 x 10 13 


1.26 x 10 14 


(i~ DAR 


v^ and v e 


5.96 x 10 10 


4.87 x 10 10 


1.08 x 10 11 


7T+ DIF 




1.37 x 10 12 


8.26 x 10 11 


2.20 x 10 12 


TT" DIF 




1.45 x 10 11 


1.11 x 10 11 


2.56 x 10 11 


7T+ DIF 




5.56 x 10 8 


5.01 x 10 s 


1.06 x 10 9 


H+ DIF 




4.13 x 10 9 


2.44 x 10 9 


6.57 x 10 9 



TABLE IV: Cross section uncertainties for the neutrino reactions with two-body final states that 
occur in LSND. The cross sections for these processes are known accurately because either related 
measurements can be used to constrain the matrix elements or only fundamental particles are 
observed. Also shown are the corresponding neutrino flux constraints. 



Process 


a Constraint 


a Uncertainty 


Flux Constraint 


ve — > ve 


Standard Model Process 


1% 


fi + — > v e v^e + DAR 


12 C(ve,e-) 12 N g . s . 


12 AT 


5% 


/i + — > v e v^e + DAR 


12 C( W -) 12 ^M. 


12 AT 


5% 


7T+ -► Vll n+ DIF 




neutron decay 


5% 


7T~ — > V^jJL - DIF 
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TABLE V: The average efficiencies for electrons in the fiducial volume with energies in the range 
20 < E e < 60 MeV. 



Criteria 


Efficiency 


Electron Reduction 


Energy > 15 MeV 


1.00 


Veto Hits < 4 


0.98 ±0.01 


No Laser Tag 


1.00 


Loose Electron PID 


0.96 ±0.01 


Vertex > 10 cm from PMTs 


1.00 


Cosmic Muon Cut 


0.92 ±0.01 


Electron Selection 


At past > 12/iS 


0.96 ±0.01 


At future > 8fJ,S 


0.99 ±0.01 


No bottom veto hit 


1.00 


-1.5 < xL < 0-5 


0.84 ± 0.01 


0-3 < xtt < 0-65 (1993 only) 


0.98 ±0.01 


85ns < t cvcnt < 210ns 


1.00 


At b v %t > 30ns 


0.97 ±0.01 


D > 35 cm 


0.88 ± 0.02 


N 7 < 1, E > 60 


1.00 


iV 7 < 2, E < 60 


1.00 


Deadtime 


DAQ & Tape Deadtime 


0.96 ± 0.02 


Veto Deadtime 


0.76 ± 0.02 


Total 


0.42 ± 0.03 
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TABLE VI: Event categories used to determine the number of events from standard neutrino 
processes. 



Category 


Past Event 


Primary Event 


Future Event 


e 








e 




f e 


12 A decay 


e 7 






n capture 


/j, e 


A* 


e(muon decay) 




fj, e {3 


A* 


e(muon decay) 


12 A decay 


e 7 P 




f e 


accidental 7 + 12 .A decay 


H e 7 


A* 


e(muon decay) 


n capture 


e no (3 
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TABLE VII: Parameters adjusted during the least squares fit procedure, along with the fitted 
correction values, central correction values, and nominal parameter values. 



Parameter 


Fitted Correction Value 


Central Correction Value 


Nominal Parameter Value 


Flux Parameters 


$DIF 


0.88 ±0.09 


1.00 ± 0.15 


0.22 x W 13 u/cm 2 


$DAR 


1.01 ±0.05 


1.00 ± 0.07 


12.6 x 10 n u/cm 2 


~t" ratio 

7T ' 


0.90 ±0.19 


1.00 ± 0.10 


0.12 


Cross Section Parameters 


a{yfC - M - 12 A*) 


0.68 ±0.23 


1.00 ± 0.25 


15.2 x 10~ 40 cm 2 




0.97 ±0.05 


1.00 ± 0.05 


4.9 x 10~ 40 cm 2 


<t(^ 12 C -> e" 12 iV 9 . s .) 


1.01 ±0.05 


1.00 ± 0.05 


9.2 x 10~ 42 cm 2 


a(i/ e 12 C -> e" 12 A*) 


1.02 ±0.13 


1.00 ± 0.25 


4.1 x 10" 42 cm 2 


a(z/ e 13 C -► e" 13 A) 


0.93 ±0.28 


1.00 ± 0.30 


0.53 x 10~ 40 cm 2 


Efficiency Parameters 


e M 


1.00 ±0.06 


1.00 ±0.07 


0.93 




1.00 ±0.04 


1.00 ±0.07 


0.65 


e e 


1.00 ±0.05 


1.00 ±0.07 


0.42 


e 7 


0.91 ±0.03 


1.00 ±0.07 


0.60 


duty ratio 


0.95 ±0.03 


1.00 ± 0.03 


0.060 
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TABLE VIII: The estimated number of events in the 20 < E e < 60 MeV energy range due to 
100% Un — > v e transmutation and to the two beam-related backgrounds with neutrons, \i~ decay 
at rest in the beam stop followed by 9 e p — > e + n scattering in the detector and ir~~ decay in flight in 
the beam stop followed by v^p — > fi + n scattering. The 7r~ DIF background includes contributions 
from v^C — > fi + nX and v^C — ► [i~nX scattering, as well as a small P e background from 7r~ and 
/i~ DIF. The events must satisfy the electron selection criteria, but no correlated 7 requirement is 
imposed. 



Neutrino Source 


Reaction 


Number of Events 


//+ DAR 
\T DAR 
7T- DIF 


100% ^ -► v e 
i> e p — > e + n 
W -»• 


33300 ± 3300 
19.5 ±3.9 
10.5 ±4.6 



TABLE IX: The correlated and accidental 7 efficiencies for different i? 7 selections. The systematic 
uncertainty of these efficiencies is estimated to be ±7% of their values. 



Selection 


Correlated 7 Efficiency 


Accidental 7 Efficiency 


i? 7 > 1 


0.51 


0.012 


R-y > 10 


0.39 


0.0026 


i? 7 > 100 


0.17 


0.0002 
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TABLE X: Numbers of beam-on events that satisfy the selection criteria for the primary — > v e 
oscillation search with TL, > 1, TL, > 10, and i? 7 > 100. Also shown are the beam-off background, 
the estimated neutrino background, the excess of events that is consistent with neutrino oscillations, 
and the probability that the excess is due to a statistical fluctuation. 



Selection 


Beam-On Events 


Beam-Off Background 


v Background 


Event Excess 


Probability 


i? 7 > 1 


205 


106.8 ±2.5 


39.2 ±3.1 


59.0 ± 14.5 ±3.1 


7.8 x 10~ 6 


i? 7 > 10 


86 


36.9 ±1.5 


16.9 ±2.3 


32.2 ±9.4 ±2.3 


1.1 x 10" 4 


i? 7 > 100 


27 


8.3 ±0.7 


5.4 ± 1.0 


13.3 ±5.2 ±1.0 


1.8 x 10" 3 



TABLE XI: The number of excess events in the 20 < E e < 60 MeV energy range, together with 
the corresponding oscillation probability if the excess is due to — > v e oscillations. Also shown 
are the results from the analysis of the 1993-1995 data sample II]. 



Analysis 


Excess Events 


Oscillation Probability 


Present Analysis (1993-1998) 
Previous Analysis (1993-1995) 


87.9 ±22.4 ±6.0 
51.01^5 ± 8.0 


(0.264 ±0.067 ±0.045)% 
(0.31 ±0.12 ±0.05)% 
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TABLE XII: The oscillation probabilities obtained with various selections. The S > 0.5 selection 



was used in the previous analysis M. The nominal values are shown in Table XI 



Selection 


Oscillation Probability 


Nominal 


(0.264 ± 0.067 ± 0.045)% 


Nominal ± Atp as j > 20//S 


(0.220 ± 0.064 ± 0.045)% 


Nominal + Veto Hits < 2 


(0.303 ± 0.074 ± 0.045)% 


Nominal + -1 5 < vl j < 


fO 304 ± 077 ± 045)% 


Nominal + D > 50 cm & Y > —50 cm 


(0.252 ±0.071 ±0.045)% 


Nominal + D > 75 cm 


(0.222 ±0.074 ±0.045)% 


Nominal + Y > -120 cm 


(0.239 ± 0.061 ± 0.045)% 


Nominal + At past > 15.2^s &Y> -120 cm 


(0.193 ±0.055 ±0.045)% 


Nominal + S > 0.5 


(0.293 ±0.069 ±0.045)% 



TABLE XIII: Number of beam on-off excess events that satisfy the selection criteria for the primary 
Vfj, — > v e oscillation search with 1 associated 7 and with > 1 associated 7. (An associated 7 is defined 
to have i? 7 > 10.) The excess of events with > 1 correlated 7 is approximately zero, which is what 
is expected for the reaction T> e p — > e + n. 



Energy Selection 


1 Associated 7 


> 1 Associated 7 


20 < E e < 60 MeV 
36 < E e < 60 MeV 


49.1 ±9.4 
28.3 ±6.6 


-2.8 ± 2.4 
-3.0 ± 1.7 
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TABLE XIV: Number of beam on-off excess events that satisfy the selection criteria for the primary 
V\i ~^ & e oscillation search with 36 < E e < 60 MeV and with > 11 "lookback" hits in the 0-3 fis 
and 3-6 (is intervals. Results are shown for events with i? 7 > and for events with i? 7 > 10. 
The number of excess events in each 3fis interval is consistent with the probability of having an 
accidental lookback in the time interval. 



i? 7 Selection 


— S/js 


3 — 6/is 


Events Expected Due to Accidentals 


i? 7 > 
i? 7 > 10 


11.5 ±6.3 
1.7 ±1.4 


7.8 ±5.9 
0.5 ± 1.0 


10.8 ±2.2 
1.6 ±0.4 



TABLE XV: The number of excess events in the 60 < E e < 200 MeV energy range, together with 
the corresponding oscillation probability if the excess is due to i/„ — > v e oscillations. Also shown 
are the results from the higher precision analysis of the 1993-1995 data sample [||. 



Analysis 


Excess Events 


Oscillation Probability 


Present Analysis (1993-1998) 
Previous Analysis (1993-1995) 


8.1 ± 12.2 ± 1.7 
18.1 ±6.6 ±4.0 


(0.10 ±0.16 ±0.04)% 
(0.26 ±0.10 ±0.05)% 
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TABLE XVI: The eight contributions to the (sin 2 29, Am 2 ) £ fit from all of the signal and 
background processes. Also shown are the fitted number of events at the best fit point of 
(sin 2 2e,Am 2 ) best ^ fu = (0.003, 1.2eV 2 ). 



C Contribution 


Signal or Background Source 


Process 


Fitted Number of Events 


1 


Vfl -> v e 


P e p — ► e + n 


89.5 


2 


BUB 




3664.6 


3 


DAR v e 


^ e 12 C -> e"iV ff . a . 


1865.0 






i/ e 12 C -> e~ A* 








i/ e 13 C -> e~A 








z^e — ► z/e 




4 


DIF Vlk 


i/„C -> fi-N* 


37.3 






v^C -> n~N g . s , 




5 


DIFF M 




5.9 


















6 


DAR V e (n~ DAR) 


V e p — ► e + n 


16.7 


7 


V(i -»• v e 


zv e C -► e~A 


6.1 


8 


DIF 7r + — > v e and fi + — > ^ e decay 


zv e C -► e~A 


11.9 



TABLE XVII: A comparison of the LSND and KARMEN experiments. 



Property- 


LSND 


KARMEN 


Proton Energy 


798 MeV 


800 MeV 


Proton Intensity 


1000 (jA 


200 fiA 


Duty Factor 


6 x 10~ 2 


1 x 10~ 5 


Total Mass 


167 t 


56 t 


Neutrino Distance 


30 m 


17.5 m 


Particle Identification 


YES 


NO 


Energy Resolution at 50 MeV 


6.6% 


1.6% 
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Beam Stop 




LSND Detector Water Plug Electronics 
and Veto System Caboose 

< >~ 

30 m 



FIG. 1: The layout of the LSND detector and the A6 beam stop area. 




Upper Perimeter of Containment 
Box Shown in Phantom 



PLAN VIEW, NEUTRINO SOURCE □ High z Material 



FIG. 2: The layout of the A6 beam stop, as it was configured for the 1993-1995 data taking. 
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FIG. 3: The decay-at-rest neutrino fluxes averaged over the detector. 
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FIG. 4: The decay-in-flight neutrino fluxes averaged over the detector. 
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FIG. 5: The electron and (5 energy distributions and the time between the electron and /?, At, 
for 12 C(v e ,e~) 12 N g , s . scattering events. 



43 



Beam Excess 
ve elastic 

^ 12 C(v e ,e) 12 N 




%400 



E e MeV 




1 -0.8 -0.6 -0.4 -0.2 



0.2 0.4 0.6 0.8 1 

cosQ,, 



FIG. 6: The energy and angular distributions for inclusive electron events. E e is the electron 
energy and 6 U is the angle between the incident neutrino and outgoing electron directions. Neutrino- 
electron elastic scattering events are clearly seen near cos9 u ~ 1. 
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FIG. 7: The angular distribution (top plot) and the energy distribution (bottom plot) for neutrino- 
electron elastic scattering events with cos9 u > 0.9 and with 12 C(v e , e~) 12 N g , s . events removed. 
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FIG. 8: The muon and (3 energy distributions (electron energy equivalent) and the time between 
the muon and (5 for 12 C(v ll ,[jL~~) l2 Ng. s . scattering events. 
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FIG. 9: The Michel electron and muon energy distributions (electron energy equivalent), the time 
between the muon and electron, At, and the distance between the reconstructed electron position 
and muon position, Ar, for v^C — > (J,~N, 9^C — > n + B, and u^p — > // + n inclusive scattering events. 
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Correlated Distribution 



Uncorrelated Distribution 




FIG. 10: Distributions for correlated 2.2 MeV 7 (solid curves) and accidental 7 (dashed curves). 
The top plot shows the distance between the reconstructed 7 position and positron position, Ar, 
the middle plot shows the time interval between the 7 and positron, At, and the bottom plot shows 
the number of hit phototubes associated with the 7, Nh its . 
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FIG. 11: The i? 7 distribution for v e C — > e~N g ^ s . exclusive events, where the N g . s . (3 decays. The 
distribution is consistent with a pure accidental 7 shape. 
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FIG. 12: The -R 7 distribution for vjC — > fi N, vJJ — > n + B, and v^p — > /i + n inclusive scattering 
events. 
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FIG. 13: The individual 7 distributions from i/„C — > fJ-^N, u^C — > ^ + B, and z^p — > 
scattering for events with i? 7 > 1 (left side) and i? 7 < 1 (right side). The top plots show the 
distance between the reconstructed 7 position and positron position, Ar, the middle plots show 
the time interval between the 7 and positron, At, and the bottom plots show the number of hit 
phototubes associated with the 7, N^its- 
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FIG. 14: The distribution for events that satisfy the selection criteria for the primary 
oscillation search. 
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FIG. 15: The individual 7 distributions for events that satisfy the selection criteria for the primary 
- * v e oscillation search with R~ > 1 (left side) and i? 7 < 1 (right side). The top plots show the 
distance between the reconstructed 7 position and positron position, Ar, the middle plots show 
the time interval between the 7 and positron, At, and the bottom plots show the number of hit 
phototubes associated with the 7, Nhits- 
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FIG. 16: The energy distribution of the 1993-1998 data sample for events with i? 7 > 10. The 
shaded region shows the expected distribution from a combination of neutrino background plus 
neutrino oscillations at low Am 2 . 
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FIG. 17: The energy distribution of the 1993-1995 (top plot) and 1996-1998 (bottom plot) data 
samples for events with R-y > 10. The shaded region shows the expected distribution from a 
combination of neutrino background plus neutrino oscillations at low Am 2 . 



55 




FIG. 18: The cos 6 U distribution for events with > 1 and 36 < E < 60 MeV. The shaded 
region shows the expected distribution from a combination of neutrino background plus neutrino 
oscillations at low Am 2 . 
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FIG. 19: The spatial distributions for events with i? 7 > 10, 20 < E e < 60 MeV, and D > 10 cm. 
The shaded region shows the expected distribution from a combination of neutrino background 
plus neutrino oscillations at low Am 2 . 
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FIG. 20: The spatial distributions for events with i? 7 > 10, 20 < E e < 60 MeV, and D > 35 cm. 
The shaded region shows the expected distribution from a combination of neutrino background 
plus neutrino oscillations at low Am 2 . 
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FIG. 21: The scatter plots of the x-y and y-z spatial distributions for events with R~ > 10, 
20 < E e < 60 MeV, and D > 35 cm. Beam-on and beam-off events are shown separately. The 
ratio of the dot area in beam-off plots, to the dot area in the beam-on plots, is equal to the duty 
ratio. This gives the appropriate scale for the beam-off subtraction. 
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FIG. 22: The particle identification, x'tot> distribution for events with i? 7 > 10, 20 < E e < 60 
MeV, and D > 35 cm. The shaded region shows the expected distribution from a combination 
of neutrino background plus neutrino oscillations at low Am 2 . Oscillation candidate events are 
required to staisfy the requirement — 1.5 < x'tot < 0.5. 



60 



3 

cb 35 
Uj 

cti 

CD 

*25 



20 



15 



10 







-0.5 0.5 1 



1.5 2 2.5 3 3.5 

N veto 



FIG. 23: The veto hit distribution for events with i? 7 > 10 and 20 < E e < 60 MeV. The data 
agree well with the distribution from u e C — > e~N g . s , scattering (shaded histogram), where the 
reaction is identified by the iV s . s . [3 decay. 
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FIG. 24: The L v /E v distribution for events with R y > 10 and 20 < E e < 60 MeV, where L v is 
the distance travelled by the neutrino in meters and E v is the neutrino energy in MeV. The data 
agree well with the expectation from neutrino background and neutrino oscillations at low Am 2 . 
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FIG. 25: This figure shows the time to the previous event for i? 7 > 10 electron events prior to 
applying the At pas t > 12^s selection. In the upper graph the beam excess events are in agreement 
with our expectations for v^p — ► /j, + n and C(z/«, fx + ) 12 B* processes. With the same data on a 
smaller vertical scale, the bottom graph shows events with accidental past activities, in agreement 
with expectations from random cosmic ray backgrounds and beam related backgrounds. Note that 
most of the oscillation candidate events have no past activity, and therefore do not appear in these 
graphs. 
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FIG. 26: Favored regions in the (sin 2 28, Am 2 ) plane at 90% CL. The Feldman-Cousins, Bayesian, 
and constant-slice methods all give about the same result. 
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FIG. 27: A (sin 2 26, Am 2 ) oscillation parameter fit for the entire data sample, 20 < E e < 200 
MeV. The fit includes primary 9^ — > v e oscillations and secondary Uu — ► u e oscillations, as well 
as all known neutrino backgrounds. The inner and outer regions correspond to 90% and 99% CL 
allowed regions, while the curves are 90% CL limits from the Bugey reactor experiment, the CCFR 
experiment at Fermilab, the NOMAD experiment at CERN, and the KARMEN experiment at 
ISIS. 
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FIG. 28: The E P , z, cos^jy, and FLy projections from the 4-dimensional (sin^ 2$, Attt,^) likelihood 
fit. The points with error bars are the data. 
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